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Abstract. We report the results of applying the 2MASS Tully-Fisher (TF) relations to study the galaxy bulk 
flows. For 1141 all-sky distributed flat RFGC galaxies we construct J,H,K S TF relations and find that Kron J/ e 
magnitudes show the smallest dispersion on the TF diagram. For the sample of 971 RFGC galaxies with V-zk < 
18000 km s~ we find a dispersion otf = 0.42 m and an amplitude of bulk flow V= 199±61 km s" 1 , directed 
towards I = 301° ± 18°, b = —2° ± 15°. Our determination of low-amplitude coherent flow is in good agreement 
with a set of recent data derived from EFAR, PSCz, SCI/SCII samples. The resultant two- dimensional smoothed 
peculiar velocity field traces well the large-scale density variations in the galaxy distributions. The regions of large 
positive peculiar velocities lie in the direction of the Great Attractor and Shapley concentration. A significant 
negative peculiar velocity is seen in the direction of Bootes and in the direction of the Local void. A small positive 
peculiar velocity (100 - 150 km s _1 ) is seen towards the Pisces-Perseus supercluster, as well as the Hercules - 
Coma - Corona Borealis supercluster regions. 
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•'-h . 1. Introduction al., 1999) with the objects from " The Extended Source 

/S ' Catalog", XSC (Cutri et al., 1998; Jarrett et al., 2000). 

H i Since the pioneering work by Aaronson, Mould, Huchra of the total number of 42 3 6 RFGC galaxies, 2996 galax- 

and collaborators (1979, 1980, 1982), the infrared Tully- ies ( e _ g _ n% ) have been detected in the J, H, if s -bands. 

Fisher relation (IRTF) has been widely used in the study The RFGC cata i ogue was crea ted using the material of 

of galactic bulk flows on different scales. A rather complete the photographic sky surveys POSS-I and ESO/SERC. In 

review on the IRTF application is given in the thesis of p aper j we analyzed in detail the 2MASS-characteristics 

Bamford (2002). of flat g a i ax i es entered in the XSC. To build the B-, I- 

The cosmic flow investigations require both enormous , H-, K s - TF relations , we used: B t - magnitudes 

and homogeneous samples. One way to increase the num- from RFGC, calculated from angular diameters, taking 

ber of observational data is the sample combination proce- into account the galaxy surface brightness, as well as type, 

dure. For example, making the Mark III catalogue, Willick and other galaxy parameters; the total /-magnitudes from 

et al. (1997) compiled the infrared magnitudes and line Mathewson & Ford (1996) and Haynes et al. (1999); the 

widths from various sources and reduced them to a uni- isophotal J-, H-, K s - magnitudes measured in elliptical 

form system. apertures at a level of K s =20 mag/arcsec 2 . For a sample 

The appearance of the complete and homogeneous of 436 flat galaxies with this set of magnitudes and known 

2MASS survey Skrutskie et al. 1997 opens up new op- radial velocities and HI line widths we obtained the slope 

portunities for cosmic flow study by IRTF. In our pre- of the linear TF regression, increasing from 4.9 in the B- 

vious work (Karachentsev et al., 2002, hereafter Paper band to 9.3 in the K s - band. The derived scatter on the 

I) we identified the spiral edge-on galaxies from "The TF diagram did not show a significant decrease from the 
Revised Flat Galaxy Catalogue" , RFGC (Karachentsev et 
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blue to infrared band, and after excluding of dwarf galax- 
ies reached ~ 0.6 mag. 

In the present work we study in dipole approximation 
the bulk motion of flat galaxies from a homogeneous sam- 
ple - the RFGC catalogue using the J-, H-, K S -TF rela- 
tions based on the 2MASS photometry data. We show that 
the scatter on the TF diagram can be diminished signifi- 
cantly by two factors: a cleaning of the initial sample and 
including additional photometric parameters in the sim- 
ple TF relation. As a result, for 971 all-sky distributed flat 
galaxies the dipole solution is: V=199 km s _1 , I = 301°, 
b = —2°. The smoothed peculiar velocity field repeats, as 
a whole, the large-scale distribution of flattened galaxies 
from the 2MASS. 

2. The sample 

2.1. The HI line widths (the sample RFGC-W 50 ). 

In Paper I we presented the radial velocities of 1772 RFGC 
galaxies. However, not all of these 1772 objects have the 
HI line width measurements. To compile the list of the 
RFGC galaxies with known estimates of W50 (measured 
directly or calculated from rotation curves), we use the 
following sources: 

a) The list of flat galaxies with known velocities and 
line widths (Karachentsev et al., 2000a). This compila- 
tion consists of several subsamples observed by Giovanclli 
et al., 1997a, Makarov et al., 1997a, b; 1999, 2001, 
Mathewson et al. (1992), Mathewson and Ford (1996), and 
Matthews and van Driel (2000). b) Our identifications of 
southern RFGC galaxies with the HIPASS survey sources 
(Karachentsev & Smirnova, 2002), and with "A Catalog of 
Hl-Sclcctcd Galaxies from the South Celestial Cap Region 
of Sky" (Kilborn et al., 2002). c) The last version of LEDA 
database (Paturel et al. 1996). d) Unpublished data on the 
HI observations of RFGC galaxies at the Effelsberg and 
Nangay radio telescopes (Huchtmeier et al., 2003). The 
sources a),b),c), and d) contain, respectively, 78%, 2%, 
11%, and 9% of the whole sample. A total of 1653 sets 
of radial velocities, Vh. and line widths, Wso, were used in 
the initial list, including some multiple observations of one 
and the same galaxy. To check the "best" estimate among 
double and triple ones, we determine in a dipole approx- 
imation the distances Hr, using the TF relation "linear 
diameter - line width" from Karachentsev et al. (2000b). 
The observed radial velocities were reduced to the cosmic 
microwave background 3K system according to Kogut et 
al. (1993), and the observed line widths were corrected for 
cosmological broadening and turbulence following Tully & 
Fouque (1985). No corrections for inclination were made 
because the RFGC galaxies with their apparent optical 
axis ratio a/6 > 7 are by definition very much inclined to 
the line of sight (i > 82°). The galaxy peculiar velocities 
were calculated as V pec = Vzk — Hr. Here and hereafter 
the inferred distance Hr is expressed in km s . We re- 
tained in the RFGC-W50 sample only the galaxies from 
multiple observations whose peculiar velocities were mini- 



mal. In all the cases we took into account the galaxy mor- 
phology, a possible confusion from a near neighbour, the 
signal/noise ratio etc. After excluding multiple measure- 
ments, 1617 RFGC galaxies with known radial velocities 
and line widths were entered in the RFGC-W50 sample. 

2.2. The 2MASS photometry ( the sample 
RFGC-2MASS) 

To compile this sample, we performed the cross- 
identification between RFGC and XSC catalogues. The 
initial file contains 3001 lines with 2MASS data. We used 
for the processing the following characteristics (Jarrett et 
al., 2000, 2003): 

r2o - major isophotal radius in arcsec, neasured at the 
20 mag/arcsec 2 level in the K s band via photometry in 
elliptical isophotes; 

r/e - fiducial elliptical Kron radius in arcseconds; 

r ext - radius of the "total" aperture in arcseconds; 

J2o/ e , #20/e, K20fe - isophotal fiducial elliptical- 
aperture magnitudes in corresponding bands measured at 
the K s - band fiducial 20 mag/arcsec 2 isophotal radius; 

Jf e , Hf e , Kf e - Kron fiducial elliptical-aperture mag- 
nitudes measured at the K s - band fiducial elliptical Kron 
radius; 

J ext , H ext , K ext — integral "total" magnitudes as de- 
rived from the isophotal magnitudes ( J20, H2o, K20) and 
the extrapolation of the fit to the radial surface bright- 
ness distribution. The extrapolation (r ext ) is carried out 
to roughly four times the disk scale length. (Details are 
given in Jarrett et al. 2003); 

Jhl - J-band "effective" surface brightness (at J- band 
half-light "effective" radius); 

Jcdex - J-band concentration index (3/4 vs. 1/4 light 
radius); 

sba - axis ratio (b/a) for the J + H + K s combined 
image ("super" coadd). 

A comparison of two lists, RFGC-2MASS and RFGC- 
W50. gives 1215 common galaxies, including 68 galaxies 
with multiple 2MASS estimates. The selection of the best 
data among the duplicate ones is made by comparing the 
deviations from the simple linear TF relation: 

M = Ci + C 2 logW c , (1) 

where 

M = m-A-25-5log(V 3K /H ), (2) 

A is extinction and Ho=75 km s~ 1 Mpc _1 . It should 
be noted that the equation (2) is known to be incorrect 
at small distances since nearby galaxies around the Local 
Group do not need to have an additional ~ 600 km s" 1 
correction subtracted from their apparent velocities to be 
placed in the proper Hubble flow. Altogether, 1141 galax- 
ies were entered in the joint RFGC-W50-2MASS sample 
after the excluding procedure. Some parameters of dif- 
ferent observables are given in Table 1, where <r is the 
standard deviation. 
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Table 1. Statistical characteristics of the flat galaxy sam- 
ple 



N=1141 


min 


max 


mean 


a 


J 20 / e , mag 


7.5 


15.9 


13.1 


1.21 


H 20 fe, mag 


6.6 


15.4 


12.3 


1.23 


K 2 ofe, mag 


6.3 


14.9 


12.0 


1.23 


Jf e , mag 


7.4 


15.2 


12.9 


1.10 


Hf e , mag 


6.6 


14.7 


12.1 


1.13 


K fe , mag 


6.3 


14.4 


11.8 


1.13 


Jext, mag 


7.4 


15.9 


12.9 


1.13 


H ext , mag 


6.5 


15.3 


12.1 


1.16 


K ext , mag 


6.2 


14.8 


11.8 


1.18 


lg(r ext , arcsec) 


0.82 


2.36 


1.57 


0.22 


lg(r/ e , arcsec) 


0.81 


2.08 


1.50 


0.21 


lg(»*20/e, arcsec) 


0.70 


2.04 


1.36 


0.25 


Jhl, mag/arcsec^ 


16.8 


20.6 


19.0 


0.68 


Jcdex 


1.89 


8.69 


3.99 


0.92 


sba 


0.10 


1.00 


0.24 


0.097 


lg(W 50 , km s" 1 ) 


1.36 


2.89 


2.50 


0.163 


Wso, km s _i 


23 


782 


335 


112 


lg(V 3K , km s" 1 ) 


2.24 


4.38 


3.73 


0.269 


V3K, km s~ x 


175 


23758 


6341 


3327 



As known, the total and the Kron magnitudes 
are something like 15 - 20 % brighter than the 
isophotal magnitudes (see figures 11, 12 and 13 in 
http://spider.ipac.caltech.edu/staff/jarrett/papers/LGA 
/LG A_fig.htm). 

3. Construction of the optimal sample and 
calculation of the dipole parameters 

For each galaxy of our sample we have a set of the isopho- 
tal, Kron, and the total magnitudes in the J, H, K s - bands. 
The check of mutual correlations between all the magni- 
tudes shows that they are well correlated (Table 2). 

As an illustration of tight mutual relationship between 
the 2MASS magnitudes, in Fig.l we give the regression 
of K20f e on i?20/e, consist with the NIR colors of disk 
galaxies (Jarrett 2000; see also Fig. 20 in Jarrett et al. 
2003). 

According to the data in Table 2, all three color bands 
seem to be of equal value. 

To select between the nine different kinds of magni- 
tudes, we built TF relation (1) for each of them and cal- 
culated the characteristics of bulk motion in a dipole ap- 
proximation. At this stage the apparent magnitudes were 
corrected for Galactic extinction as 

J c = J - 0.207 ■ A B , 

H C = H -0.132- A B , (3) 

K c = K- 0.084 -A B , 

where As is an extinction in the S-band according to 
Schlegel et al. (1988). Here we did not correct the magni- 
tudes for internal galaxy extinction, because the internal 
extinction depends on the galaxy inclination, as well as 
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Fig. 1. The relationship between i?20/e and i^20/e- 

the galaxy luminosity (Vcrhcijcn, 2001), which is a priori 
unknown. 

For each of magnitudes we use constructed TF rela- 
tions to calculate the characteristics of bulk motion in a 
dipole approximation. We adopt the simple linear model: 
Vp e c — Ve r + 5, where V - the bulk velocity to the apex, 
e r - an unit radial vector of galaxy direction. The value of 
T,S 2 with summing on the whole sample have been mini- 
mized by a least squares method. 

The results of calculations are given in Table 3. 

The columns of Table 3 denote: 

(1) - 2MASS magnitude, by which the absolute mag- 
nitude in the TF relation was calculated; 

(2) - dispersion on the TF diagram in mag; 

(3) - slope C2 in formula (1) and its statistical signifi- 
cance according to the Fisher criterion in parentheses; 

(4) - dispersion of peculiar velocities in km s , which 
includes the error in the measured distances; 

(5) - modulus of the bulk motion velocity and its error 
in km s , 

(6) , (7) - galactic longitude and latitude of the apex 
and their errors in degrees; 

(8) - significance of the vectorial dipole solution ac- 
cording to the Fisher criterion. (Note that for confidence 
probability 95%, the quantile of the Fisher distribution is 
equal to 2.6 for three degrees of freedom for numerator 
and infinity ones for denominator). 

The errors in V, I, and b were calculated by first 
calculating the diagonal elements Byv,Bu and B^ 
of covariance matrix B in the frame {ey, e;, ef,},and 
then AV = {B VV ) 1/2 ^l = arctan{(B„) 1 / 2 /V A }, Al = 
arctan{( J B bb ) 1 /2/T/}. 

As seen from Table 3, for every photometric magni- 
tude, the dispersion in the TF fit is too large, nearly 1.5 
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Table 2. Some parameters of mutually correlated linear dependence y = kx + c between different magnitudes for the 
sample of 1141 galaxies 



X 


V 


(J fwd 




(J orth 


r 


k 


c 






0.119 


0.117 


0.083 


0.9954 


1.019±0.003 


-1.042±0.039 


J^Ofe 




0.162 


0.159 


0.114 


0.9913 


1.019±0.004 


-1.404±0.053 


Hon f P 




0.114 


0.114 


0.080 


0.9958 


1.000±0.003 


-0.361±0.034 


J f e 
■j j e 


f e 


0.131 


0.127 


0.091 


0.9933 


1.026±0.004 


-1.107±0.046 


Jfp 
J e 


J C 


0.173 


0.168 


0.121 


0.9882 


1.026±0.005 


-1.460±0.061 


j r - 


K fe 


0.131 


0.131 


0.093 


0.9932 


1.000±0.003 


-0.352±0.042 






0.197 


0.193 


0.138 


0.9854 


1.025±0.005 


-1.091±0.068 


jlll 


K ex t 


0.223 


0.214 


0.155 


0.9820 


1.043±0.006 


-1.648±0.077 


Hext 


K ext 


0.190 


0.187 


0.133 


0.9870 


1.018±0.005 


-0.535±0.060 


J20fe 


Jfe 


0.146 


0.161 


0.108 


0.9912 


0.905±0.004 


-1.029±0.047 


H20fe 


Hfe 


0.151 


0.165 


0.111 


0.9910 


0.911±0.004 


+0.898±0.045 


K20fe 


Kfe 


0.143 


0.157 


0.106 


0.9919 


0.911±0.003 


+0.874±0.042 


Jfe 


J ext 


0.165 


0.160 


0.115 


0.9894 


1.030±0.004 


-0.414±0.058 


Hf e 


H ex t 


0.177 


0.172 


0.124 


0.9883 


1.030±0.005 


-0.373±0.057 


Kfe 


Kext 


0.150 


0.144 


0.104 


0.9919 


1.048±0.004 


-0.544±0.047 



Notes to table: of w d, o~ inv , a ort h are the standard deviations (in mag) from the direct, 
inverse and orthogonal TF regressions, respectively; r is the coefficient of correlation; 
k and c are the coefficients of the orthogonal regression. 



Table 3. Parameters of the TF relation (1) and the dipole characteristics for 1141 flat galaxies. 





0~TF 


c 2 


Ov 


V 


/ 


b 


F 


Jfe 


0.855 


-6.24±0.15 (1653) 


1940 


416±106 


284±15 


-21±12 


5.2 


J20fe 


0.940 


-6.81±0.17 (1630) 


2177 


568±118 


292±12 


-24±10 


7.8 


J ext 


0.911 


-6.50±0.16 (1581) 


2076 


437±113 


287±15 


-17±12 


5.0 


Hfe 


0.887 


-6.50±0.16 (1668) 


2033 


448±111 


287±15 


-23±12 


5.5 


H20fe 


0.972 


-7.05±0.17 (1635) 


2267 


565±122 


293±13 


-24±10 


7.2 


Hext 


0.947 


-6.78±0.17 (1592) 


2183 


443±119 


288±16 


-24±13 


4.7 


Kfe 


0.907 


-6.68±0.16 (1686) 


2062 


340±113 


285±19 


-22±15 


3.1 


K 2 0fe 


0.990 


-7.21±0.18 (1647) 


2295 


475±124 


292±15 


~23±12 


4.9 


Kext 


0.966 


-6.99±0.17 (1625) 


2226 


378±121 


290±19 


-22±15 


3.3 



times as high as that presented in Paper I. We suggest 
therefore that the dipole parameters in Table 3 be con- 
sidered only preliminary. Due to a larger dispersion, the 
slopes of the TF relation are flatter than those obtained 
in Paper I. 

It is also seen that the Kron magnitudes (especially 
Jf e ) have a slightly lower dispersion in the TF diagram 
than isophotal as well as extended magnitudes. We built 
the diagrams of the residuals "isophotal minus Kron mag- 
nitude" depending on the absolute isophotal magnitude. 
For all three color bands in a wide range of absolute magni- 
tudes {-16 m , -20" 1 } the Kron magnitudes are on average 
~0.4 m brighter, thus recovering a significant fraction of 
the disk light that is lost in the background noise (see also 
Fig. 12 in Jarrett et al., 2003) . This can reduce partially 
the known non-linearity of the TF relation at the lower 
luminous end. 

We adopt Jf e as the basic apparent magnitude and ac- 
complish the procedure of trimming of combined RFGC- 
W50-2MASS sample. The galaxies were omitted from the 
sample by such criteria: 



1. W 50 < 100 km s- 1 (N=13). These are bluish dwarf 
galaxies with large uncertainties in both the profile width 
and the IR magnitudes. 

2. The deviation from the TF relation is more than 3cr 
(N=63). 

3. \V pec \ > 3000 km s" 1 (N=119). We assume that 
such large peculiar velocity values are due to non-physical 
reasons. 

We also eliminate 6 galaxies with V^k > 18000 km 
s _1 and (symmetrically) 10 galaxies with Hr > 18000 km 
s _1 to diminish the incompleteness of the sample at large 
distances. The process of elimination was converged after 
5 steps; altogether, taking into account the cross-sections, 
we excluded 170 galaxies. 

Excluding the dwarf galaxies with smallest veloc- 
ity widths means simultaneously eliminating the nearest 
galaxies. For new sample of 971 galaxies the minimal value 
of Vzk increases to 494 km s _1 in comparison with 175 km 
s _1 for non-cleaned sample. Note that from the total num- 
ber of 971 galaxies, about 87% galaxies have V-$k > 3000 
km s _1 , i.e. lie outside the Local supercluster. Thus, we 
believe the contribution of nearby galaxies in the bulk 
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motion parameters (when the CMB frame is used) is not 
significant. 




ii 4 ■ u to ip h ii ii ii ii 

Fig. 2. The Tully-Fishcr relation for 971 flat edge-on 
galaxies. 

The TF relation for a new sample of 971 galaxies 
as M(Jf e ) vs. log W c is given in Fig. 2. It has a slope 
Ci = —8.1 and a dispersion cttf=0.430. The bulk motion 
parameters obtained from this relation are: 

cry = 1062fcm • s~\ V = 323km • s _1 , 

I = 293°, b= -15°. 

As is seen in Fig. 2, eliminating 15% of the initial sam- 
ple using the " clean" criteria given above, greatly reduces 
the dispersion in the TF relation. 

Next, for the sample of N=971 we consider seven- 
parametric generalized TF relation, in which photometric 
characteristics and color index are entered: 

M(J fe ) = d+C 2 - lg(W c ) + C 3 ■ jhl + C 4 ■ jcdex+ 
+C 5 • ( J fe - K fe ) + C 6 ■ sba + C 7 ■ lg{r ext ) (4) 

For the seven- parametric TF relation (4) the solutions 
are: the slope C 2 = -7.6, a T F = 0.423, oy= 1044 km s -1 , 
V = 187 km s-\ I = 300°, b = -1°. The Jhl, (J /e - K fe ) 
and sba terms are statistically insignificant. In this case, 
only four terms remain in the new generalized TF relation: 

M(J fe ) = C x +C 2 ■ \g{W c ) + C 3 ■ jcdex + C 4 ■ \g(r ext ). (5) 

The new parameters are practically the same as for 
(4): zero-point G x = -1.6, slope C 2 = -7.6, C 3 = -0.046, 



C 4 = -0.44, ct T f= 0.422, o-y = 1045 km s" 1 , V = 199 km 
s" 1 , I = 301°, b = -2°. 

The significance of the jcdex term (^=4.5) means the 
TF relations depend on the Hubble type, namely the cen- 
tral bulge to disk distribution. Statistical insignificance of 
the sba term (^=0.3) indicates that internal extinction in 
the infrared does not correlate, surprisingly, with the in- 
frared axial ratio. The influence of the surface brightness 
term, jhl, in the regression (4) falls, but its role is par- 
tially evident via the \g{r ext ) term (F=17.5). In Fig. 3 we 
give distances Hr derived from TF relation (5) depending 
on radial velocities Vzk for the sample N=971. 




— 




ha m.u 



Fig. 3. The relationship between the derived distances 
Hr and radial velocities V^k- 

We also built analogous (seven- and four- parametric) 
TF relations based on the other eight 2MASS photomet- 
ric measures. They demonstrate the same tendencies: the 
terms Jhl, (Jf e — Kf e ), sba also are insignificant and the 
dipole bulk motion parameters are near to ones mentioned 
above. 

The results of calculations with nine TF relations (5) 
for the sample of galaxies with V^k < 18000 km s _1 and 
Hr < 18000 km s" 1 are collected in Table 4. The desig- 
nations are the same as in Table 3. 

We constructed also the second sample with the same 
restrictions as for the sample iV=971, but with V^k < 
12000 km s -1 and Hr <12000 km s _1 , to diminish the in- 
fluence of data incompleteness at large distances. In Table 
5 we show the same data as in Table 4 but for the sample 
of galaxies with V 3K < 12000 km s" 1 and Hr < 12000 km 
s _1 .The designations are the same as in Table 3. The elim- 
ination of 50 remote galaxies does not change the dipole 
parameters much as well as the TF parameters Ci. 

Comparing the data in Tables 3-5 allows us to con- 
clude: 
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Table 4. Parameters of TF regressions (5) for the 2MASS magnitudes and the dipole characteristics for 971 flat 
galaxies. 







c 2 


(TV 


V 


I 


b 


F 


t) to 


0.422 


-7.64 ± 0.13 (3558) 


1045 


199±61 


301±18 


-2±15 


3.5 


J20fe 


0.457 


-8.21±0.14 (3505) 


1144 


248 ± 67 


306 ± 16 


-5 ±13 


4.6 


Jext 


0.464 


-7.85±0.14 (3111) 


1153 


215 ± 67 


302 ± 18 


+7 ±15 


3.5 


Hf e 


0.437 


-7.86±0.13 (3509) 


1078 


209±64 


294±18 


-9 ±16 


3.6 


H20fe 


0.471 


-8.45±0.14 (3490) 


1167 


239±69 


304±17 


-7±16 


4.0 


H ex t 


0.485 


-8.15±0.15 (3067) 


1194 


198±70 


288±20 


0±17 


2.7 


Kfe 


0.445 


-8.14±0.13 (3635) 


1116 


153 ± 66 


294±24 


-1±20 


1.8 


K20fe 


0.483 


-8.72 ±0.15 (3552) 


1209 


201±71 


300±20 


-3±18 


2.7 


K ex t 


0.480 


-8.34 ±0.15 (3289) 


1204 


138±69 


301±28 


+ 12±24 


1.4 



Table 5. Parameters of TF regressions (5) for the 2MASS magnitudes and the dipole characteristics for 921 flat 
galaxies. 





(JTF 


c 2 


(TV 


V 


I 


b 


F 


Jfe 


0.432 


-7.64±0.14 (3033) 


1023 


226±62 


295 ± 16 


-2±13 


4.5 




0.466 


-8.20 ±0.15 (2993) 


1111 


271±67 


301±15 


-5 ±12 


5.4 


Jext 


0.473 


-7.83 ±0.15 (2657) 


1113 


254±66 


299±16 


±6 ±13 


5.0 


Hfe 


0.446 


-7.84 ±0.14 (2995) 


1045 


242 ± 64 


289 ± 15 


-12 ±14 


4.7 


H20fe 


0.481 


-8.42±0.15 (2978) 


1129 


263 ± 69 


298 ± 15 


-9 ±14 


4.9 


H ex t 


0.495 


-8.11±0.16 (2601) 


1156 


243 ± 71 


289 ± 17 


-4 ±14 


4.0 


Kfe 


0.453 


-8.10±0.15 (3101) 


1067 


195±65 


288±19 


-4±16 


3.0 


K20fe 


0.491 


-8.68±0.16 (3033) 


1151 


240 ± 70 


295±17 


-5 ±15 


3.9 


Kext 


0.487 


-8.28±0.16 (2802) 


1141 


183 ± 68 


295±21 


±6 ±18 


2.5 



1) The cleaning of the sample diminishes significantly 
the dispersion gtf (~ 2 times) and the bulk motion ve- 
locity modulus (about 1.5 - 2.5 times). 

2) In all procedures (the cleaning of the sample, the 
use of regression with other magnitudes) the galactic lon- 
gitude of apex changes little, within a 20° range. 

3) The cleaning of the sample moves the apex to the 
galactic equator. 

4. Peculiar velocity field 

Basing on four-parametric regression (5) and derived pe- 
culiar velocities, we built the peculiar velocity field of 
the RFGC galaxies for the sample 7V=921. In Fig. 4 we 
present a smoothed peculiar velocity field for them in su- 
pergalactic coordinates SGL, SGB. The solid line indi- 
cates the "Zone of Avoidance" within ±10 degrees around 
the Galactic equator, where effects from reddening, stel- 
lar contamination and incompleteness arc significant. The 
galaxies with V pec > are designated by open circles, and 
the galaxies with V pec < with filled circles. The smooth- 
ing was done with a Gaussian filter of 20° window. The 
color scale at the bottom Fig. 4 shows a range of average 
peculiar velocities from —200 km s _1 to ±500 km s _1 . 

In Fig. 5 we give the sky distribution in SGL, SGB 
of 22361 galaxies from the 2MASS survey selected by in- 
frared axial ratio criterion sba < 0.3. Bottom right the 
positions of the Great Attractor (GA), the Pisces- Perseus 
supercluster (PP), the Shapley concentration (Sh), as well 
as the Local Void (LV) and the Bootes Void (BV) are 
marked. Comparison of Figs. 4 and 5 shows that the re- 



gions with positive and negative smoothed velocities re- 
veal, as a whole, the galaxy overdensed and underdenscd 
regions. The regions with maximal positive V pec lie in the 
direction of the GA and the Shapely Concentration. We 
should draw a cautionary note that the obscuring effects 
on the Milky Way are also at its maximum towards these 
regions (e.g., Abell 3627 of the GA region is well within the 
" Zone of Avoidance" ) , rendering large regions of the sky 
incomplete in the 2MASS XSC and the RFGC. Significant 
negative V pec are seen towards the voids, and low positive 
velocities (~100 - 150 km s^ 1 ) are seen in the PP and 
Hercules - Coma - Corona Borealis superclusters. 

5. Comparison with other data 

To estimate the amount and quality of the observational 
material by which one was used to determine the bulk 
motion parameters, we introduce the value G, meaning 
the sample goodness 

G= (A7100) 1/2 At T f. (6) 

Here N is the number of galaxies under study and <jtf 
is their dispersion on the TF diagram. From the papers 
of different authors, we took the sample of spiral galax- 
ies lying mainly in sparse regions. In columns of Table 6 
we give: (1) - designation of the sample, (2) - color band 
used, (3) - number of galaxies in the sample, (4) - the 
dispersion on the TF diagram, (5) - the sample goodness, 
(6) - source of data. The upper four lines note the sam- 
ples included in the Mark III catalog with their respective 
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Fig. 4. A smoothed peculiar velocity field for 921 RFGC galaxies. The galaxies with V p , 
designated as open and filled circles, respectively. The "Zone of Avoidance" is indicated by a solid line 



> and V pec < are 



name. The data for them are taken from the Tables 1 and 
3 in Willick et al. (1997). The next three samples are inter- 
esting because of their relativelly large depth (cz ~ 10000 
km s _1 ). The sample of Bamford (2002) was selected from 
the SCI sample of Giovanelli et al. (1997b). From 782 S 
galaxies placed in the clusters or groups Bamford selected 
153 galaxies with available 2MASS photometry. Besides 
Bamford's data, we have not yet encountered published 
TF results using 2MASS. 

In last line we give the results obtained in this paper. 
As seen, the data goodness varies over a rather wide range, 
and the enormous samples have, naturally, the higher 
value of G. 

Note, that different authors use various manners for 
creating their "pure" samples. In case of our sample, we 
have not yet made a detailed analysis of the observable er- 
rors, nor reduced the HI line measurements to a common 
system, nor analyzed the Malmquist bias. We excluded the 
outlying galaxies only statistically because of their greater 
than 3cr deviations from the TF fit and large calculated 



V pec . Nevertheless, our data are of high goodness (G=7.4 
in comparison with the median 4.8 for the samples in Table 
6). As the minimal scatter on the TF diagram has a finite 
intrinsic value (0.25 - 0.30 mag), the improvement of the 
data goodness is undoubtedly associated with the signifi- 
cantly larger samples. 

Let us compare briefly our results with the recent liter- 
ature data. The most representative sample in Table 6, the 
Mark III catalog, yields the bulk motion V = 370 ± 110 
km s _1 towards I = 305°, b = +14° in the volume of 
Hr < 5000 km s" 1 (Dekel et al.,1999). Our dipole solu- 
tion is consistent with the Mark III data. Zaroubi (2002) 
presents a review of the latest results of galaxy bulk mo- 
tion measurements accomplished by different authors on 
various observational data. We compare our result with 
other dipole determinations using Fig.l and Table 1 from 
Zaroubi (2002). Our RFGC sample can be considered as 
relatively distant (its depth reaches 18000 km s _1 ). At 
the comparable distances (~100 - 150 Mpc h^ 1 ) accord- 
ing to Zaroubi compilation, the bulk motion parameters 
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Fig. 5. The sky distribution in SGL, SGB of 22361 galaxies from 2MASS survey with an axial ratio sba < 0.3. Right 
bottom the positions of the Great Attractor (GA), the Pisces-Perseus (PP), the Shapley concentration (Sh), as well 
as the Local Void (LV), the Bootes Void (BV) are marked. The "Zone of Avoidance" within Galactic latitude of ± 10 
degrees is shown by grey color. 

Table 6. Values of G for several samples of spiral galaxies 



Sample 


Band 


N 


OTF 


G 


Reference 


A82 field sample 


#-0.5 


359 


0.47 


4.0 


Willick et al.,1997 (Mark III) 


MAT field sample 


/ 


1355 


0.43 


8.6 


Mark III 


W91PP field sample 


r 


326 


0.38 


5.0 


Mark III 


CF field sample 


r 


321 


0.38 


4.7 


Mark III 


Great Wall spirals 


I 


172 


0.32 


4.2 


Dell' Antonio et al, 1996 


Spirals in clusters 


I 


522 


0.38 


6.0 


Dale et al, 1999 


Spirals in clusters 


J, H, K s 


153 


0.48 


2.6 


Bamford, 2002 


RFGC 


J, H, K s 


971 


0.42 


7.4 


This work 



were determined for the samples of distant elliptical galax- 
ies, EFAR (Colless et al., 2001), distant Abell clusters, 
LP10 (Willick, 1999), point sources from IRAS z-catalog 
(PSCz) (Saunders et al. 2000, Branchini et al. 2000), spiral 
galaxies in clusters and superclusters (Dale et al., 1999), 
and distant clusters, LP (Lauer & Postman, 1994). 

The low amplitude of bulk motion, V = 199 ± 62 km 
s _1 , determined in the present paper for 971 RFGC galax- 



ies with Hr <18000 km s^ 1 and V ZK < 18000 km s" 1 , is in 
agreement with the results V ~(0 - 200) km s _1 obtained 
for the samples EFAR, PSCz, SCI/SCII, as well as SNIa 
(Riess et al., 1997) which is consistent with the assump- 
tion of the flow field convergence to the CMB rest-frame 
at 100 h" 1 Mpc. However, our result differs strongly from 
the bulk velocity obtained from LP10 and LP samples 
(V ^700 km s _1 ). The low bulk velocity at large scales 
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has been predicted by popular theories of structure for- 
mation in the cosmological model with cold dark matter 
and cosmological constant. 

The bulk velocity apex obtained in this paper, I = 
301° ±18°, b = —2° ±15°, is located near the apex position 
(I = 282°, b = -8°) derived from SNIa by Riess et al. 
(1997) and, within the errors, lies near GA and the massive 
cluster Abell 3627. 

6. Conclusion 

We used the 2MASS TF relations to obtain parameters of 
bulk motion for flat edge-on galaxies from RFGC. Because 
about 71% of the all RFGC sample have J, H, K s mag- 
nitudes from 2MASS, the size of our sample under study 
is restricted only by the available HI line width data. At 
present, in there is a sample of 1141 all-sky distributed 
RFGC galaxies in our disposal with both 2MASS magni- 
tudes and velocity/line- width estimates. 

After excluding about 15% of the sample (dwarf galax- 
ies, very distant ones, and also the objects which have 
deviation over 3cr on the TF diagram), we built a set of 
multi- parametric TF relations using the Kron, isophotal, 
and extended J, H, K s magnitudes. The minimal disper- 
sion on the TF diagram is shown to be that for Kron Jf e 
magnitude {a TF = 0.422 m ). 

For the RFGC sample of A=971 galaxies with V 3 k < 
18000 km s _1 , the bulk velocity and apex position are: 
V = 199± 61 km s" 1 , I = 301° ± 18°, b = -2° ± 15°. 
The parameters of bulk motion change insignificantly with 
the use of other 2MASS magnitudes J20/C 
Hf e , Kf e , J ext , H ext , K ext (V changes within ± 50 km 
s _1 , I and b within ±15°). This result remains robust to 
decreasing the sample depth. In particular, for 921 RFGC 
galaxies with V 3K < 12000 km/s we obtained V = 226±62 
km s" 1 , I = 295° ± 16°, b = -2° ± 13°. Within the errors, 
our estimates of the bulk velocity are in agreement with 
the data for the EFAR, PSCz, SCI/SCII samples. Besides, 
the apex position and V are consistent with the results 
obtained for SNIa. 

The two-dimensional smoothed peculiar velocity field 
well traces the large-scale density variations in the 
galaxy distribution, e.g. Great Attractor, Pisces-Perseus 
an Hercules-Coma-Corona Borealis superclusters, Shapley 
concentration, Bootes and Local voids. 

Thus, we show by the example of the RFGC catalogue 
that the 2MASS Tully-Fisher relation can be used suc- 
cessfully in studying galaxy cosmic flows. The next obvi- 
ous step is to complete the observations of radial velocities 
and HI line widths for the remaining RFGC galaxies. 
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